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Abstract 

We examine several higher spin modes of the Poincare gauge theory (PGT) of gravity 
using the Hamiltonian analysis. The appearance of certain undesirable effects due to 
non-linear constraints in the Hamiltonian analysis are used as a test. We find that the 
phenomena of field activation and constraint bifurcation both exist in the pure spin 1 and 
the pure spin 2 modes. The coupled spin-0“ and spin-2“ modes also fail our test due to 
the appearance of constraint bifurcaHon. The “promising” case in the linearized theory 
of PGT given by Kuhfuss and Nitschcl likewise does not pass. From this analysis of these 
specific PGT modes we conclude that an examination of such nonlinear constraint effects 
shows great promise as a strong test for this and other alternate theories of gravity. 


1 Introduction 

Many alternative theories of gravity have been proposed to replace general relativity (GR); ex¬ 
perimental and observational results are applied to examine the predictions of each alternative 
and to eliminate the unsatisfactory onesl. Certain gauge theories of gravity, with their many 
ranges of parameters, not only pass all the experimental and observational tests but also agree 
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with GR to Post-Newtonian ordeiH. Consequently no discriminating experimental and obser¬ 
vational test can presently distinguish them from GR. Thus they are observationally equally 
viable. Spontaneously many people paid attention to another criteria — theoretical tests. 
Essentially, besides self-consistency, theoretical tests judge if a theory obeys commonly ac¬ 
cepted physical requirements, such as “no tachyons” (faster than light signals) and “no ghosts” 
(negative energy waves). Moreover a viable theory should have an appropriate mathematical 
structure. This includes a well-posed initial value problem: the basic requirement is the Cauchy- 
Kovalevska theorem; beyond that, the propagation of dynamic modes should be described by 
hyperbolic evolution equations with well-behaved characteristics. 

Although the metric-affine gauge theory of gravity (MAG)§ is more general (it has curva¬ 
ture, torsion and non-metricity), we here apply our Hamiltonian based analysis mainly to the 
Poincare gauge theory (PGT), an important sub-theory (with only curvature and torsion, i.e., 
vanishing nonmetricity) of the MAG. One reason for this choice is physical. In the foreseeable 
future the effect of nonmetricity is thought to be extremely hard to measure as it is expected 
to only survive at the Planck energy scale. On the other hand it is believed that the effect of 
torsion has the possibility to be detected from some strong astronomical sources. However our 
primary motivation for restricting our considerations to the PGT is that it is simpler. Before 
examining nonlinear effects in the MAG, we really need more experience. Also there have been 
many more studies of the behavior of the PGT. These works provide a support for our deeper 
investigation. 

Pioneered by Sezgin and Nieuwenhuizenll, people examined the PGT based on the “spin- 
projection operator formalism”i and the weak field approximation with positive energy and 
mass0. They put restrictions on the ranges of the PGT parameters and gave tables of “vi¬ 
able” PGT parameter ranges with well-behaved propagating modes. The positive energy test 
(PET), which excludes the cases permitting a non-trivial solution that has non-positive en¬ 
ergy, was proposed as an effective theoretical test for gravitational theories and applied to the 
PGTii. Parameter restrictions leading to the interesting purely curvature or torsion modes 
were identified^. On the other hand, Gauchy-Kovalevska and hyperbolicity conditions were 
imposed on the PGT to study the initial value problemjni; a shock wave analysis was applied 
to find the characteristics of propagating modesld. However, it should be noted that nonlinear 
effects were essentially not probed in most of these investigations. Almost all of the above 
examinations were directly or effectively working on the linearization of the PGT, except for 
the PET. 

The PET has been applied to the PGT and gave some good results with only simple ansatz 
which are used to find “bad” solution^. The nice feature is that the PET considers the full 
nonlinear theory and thus is not subject to such doubts that a result of a test derived from the 
linearized theory suffers. However, the test can be applied only if the total energy is well defined, 
hence it applies only to the asymptotically fiat or constant negative curvature configurations. 
Actually, in practice the PET has not given many new results beyond those found from the 
linearized theory tests, simply because it is not so easy to apply in general. This limitation 
provided further motivation for our nonlinear Hamiltonian investigation. 

An early warning of the difficulties (due to nonlinearity) was raised W Kopczynski0. He 
noted that the one parameter teleparallel theory (OPTP, a.k.a. NGR)Ej had predictability 
problems. Further analysis^ revealed that the problem was not generic, it occurred onlv for a 
special class of solutions. Gheng et al using a Hamiltonian approach verified this propertyEJ and 
identified the troubles with the phenomenon of constraint bifurcation: the chain of constraints 
could “bifurcate” depending on the values of the fields, so that the number and/or type of 
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constraints would depend on the values of the phase space variables. Since the OPTP is a 
sub-theory (with curvature vanishing) of the PGT, it was expected that the PGT as well as 
the MAG and other alternate geometric gravity theories were also vulnerable to this difficulty. 
Besides constraint bifurcation, it is likely that the gauge theories of gravity will encounter 
the nonlinear phenomenon of turning constraints into equations of motion as in the analysis 
by Velo and Zwanzingei@ of higher-spin helds. Also within Einstein’s theory of GR with 
sources included, the inherent non-linearity of the gravitational coupling to sources has long 
been known to have peculiar effects, including the activation of non-dynamic source degrees 
of freedom (see, in particular. Ref. 19 and §5.2 in Ref. 20). These interesting considerations 
motivated our investigation which is reported in this paper. 

Our main tool for studying these nonlinear effects is the Hamiltonian analysis using the Dirac 
constraint algorithm in the full nonlinear PGT. The Hamiltonian analysis reveals whether a 
theory is consistent. It systematically determines the constraints, the gauges, and the degrees 
of freedom for each consistent theory. The analysis also prescribes methods for reducing the 
phase space of the theory to one containing only the degrees of freedomil. Moreover it is the 
most straightforward way to hnd the apparatus of the initial value scheme for a given theory 
— the proper data, the constraints, and the evolution equations. 

We base our test of the viability of a particular case of the PGT (with certain parame¬ 
ter choices) by checking if the whole structure of the theory stays the same before and after 
linearization. The change of the structure of the theory can be measured by studying the non¬ 
linearity in the Poisson matrix of the constraints. In practice this was not so complicated, for 
in all the cases we have examined it was sufficient to consider the Poisson matrix of the primary 
constraints. In Ref. 21 we analyzed two modes of PGT with propagating spin-0 torsion. In 
this paper we continue our examination of the PGT considering now the simple spin-1 modes, 
simple spin -2 modes, coupled spin- 0 “ and spin- 2 “ modes, and the modes of the preferred pa¬ 
rameter choice identihed by Kuhfuss and NitschB. We hnd that two phenomena, i.e., constraint 
bifurcation and field activation, could happen in all of these caseJ^. These phenomena prevent 
these cases from being viable. 

This paper is organized as follows. In section 2 we review the basic elements of the PGT 
and introduce its Lagrangian and held equations. In section 3 we use the Dirac theory for 
constrained Hamiltonian systems in the neat “if” constraint formulation developed by Blago- 
jevic and Nikolidilil. The primary constraints, including ten “sure” primary constraints and 
thirty so-called “primary if-constraints”, are given. The total Hamiltonian density, including 
the canonical Hamiltonian density and all possible primary constraints, is presented. We sum¬ 
marize the results and the viable parameter combinations of the linearized PGT in section 4. 
The viable parameter combinations in this section are a guide to search for the viable modes of 
the full PGT. In section 5 we study the viability of the PGT of the simple spin-1 modes, simple 
spin -2 modes, coupled spin- 0 “ and spin- 2 “ modes, and give an argument to explain the failure 
of the parameter choice by Kuhfuss and Nitsch. In the hnal section we present our conclusions. 

Throughout the paper our conventions are basically the same as Hehl’s in Ref. 25. We have 
made a few adjustments to accommodate the translation of the Hamiltonian “if” constraint 
formalism to these conventions. The Latin indices are coordinate (holonomic) indices, whereas 
the Greek indices are orthonormal frame (an-holonomic) indices. The hrst letters of both 
alphabets (a, 6 , c,...; a,/?, 7 ,...) run over 1, 2, 3, whereas the later ones run over 0, 1, 2, 
3. Furthermore, =diag(—, -|-, -|-, -[-); is the completely antisymmetric tensor with 

g0i23 _ rjpg meaning of a bar over a Greek index is adopted from Ref. 24. 
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2 Poincare Gauge Theory of Gravitation 

In the PGT there are two sets of gauge potentials, the orthonormal frame held (tetrads) 
and the metric-compatible connection , which are associated with the translation and the 
Lorentz subgroups of the Poincare gauge group, respectively. The associated held strengths are 
the torsion 

Tif = + ri.i/ei,/), (1) 

and the curvature 

which satisfy the Bianchi identities 

'^[iTjkf = Riijkf-, (3) 

= 0. (4) 


The conventional form of the action, which is invariant under the Poincare gauge group, has 
the form 

S = j (fxe{LM + Tg), (5) 

where Lm stands for the matter Lagrangian density (which determines the energy-momentum 
and spin source currents) , Lg denotes the gravitational Lagrangian density, and e = det{ei^). 
In this paper we are concerned with the gravitational propagating modes, hence we omit the 
matter Lagrangian density, so Lq is considered as the source-free total Lagrangian. Varying 
with respect to the potentials then gives the (vacuum) field equations. 


with the held momenta 


- e; = 0 , 


deLp _ ^ deLg 
ddjCF 

deLp _ deLp 


and 





e\eLp - 


( 6 ) 

(7) 

( 8 ) 
(9) 


( 10 ) 

( 11 ) 


The Lagrangian is chosen (as usual) to be at most of quadratic order in the held strengths, 
then the held momenta are linear in the held strengths: 




/if/ 



( 12 ) 


(13) 
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ik).. (fc).. 

the three and the six are the algebraically irreducible parts of the torsion and the 

curvature, respectively. The reciprocal frames e\ and satisfy e\ei‘' = 5^'^ and = 

h/; the coordinate metric is defined by gij = eik’ej'^ri^y. I and k are constants; the and 

bk are free coupling parameters. Due to the Bach-Lanczos identity only five of the six b^s 
are independent. Oq is the coupling parameter of the scalar curvature R := For the 

Hamiltonian formulation we associate the canonical momenta with certain components of the 
covariant field momenta: 


TT 


- TJ- ZO 

- 7 


TT 


fii^ 


- IT ZO 

— • 


(14) 


3 Total Hamiltonian and Primary If-constrains 

The Hamiltonian analysis is our main tool for studying the PGT. We present here the primary 
constraints and the total Hamiltonian density in terms of the decomposition of the canoni¬ 
cal variables. As indicated in Table |l|, in addition to certain sure constraints associated with 
the basic local Poincare gauge symmetry of spacetime, for each of the kinetic critical param¬ 
eter combinations which vanishes there is an extra degeneracy. Such extra degeneracies are 
associated with extra primary constraints which can lead to secondaries etc. Essentially we 
follow the techniques developed by Blagojevic and NikoliCjWho investigated systematically all 
of these possibilities of the critical parameter combinations^cll. They adopted the if-constraint 
technique which avoids getting stuck in unnecessary detail and still identihes all possible con¬ 
straints. 

The sure primary constraints are 

A ^ 0, ^ 0. (15) 

These constraints reflect the fact that the torsion and the curvature are dehned as the anti¬ 
symmetric derivatives of and they do not involve the “velocities” and One 

will obtain further primary “if”-constraints if the Lagrangian density is singular with respect 
to any parts of the remaining “velocities”, and Ra^'^■ The total Hamiltonian density is 
formally of the form 

Rt = Kan + (16) 

where the (pA are the primary “if” constraints and the u’s denote the associated Lagrange 
multipliers] 'Hcan stands for the canonical Hamiltonian density which will be specified below. 
The primary constraints (^) are first-class, i.e., and are unphysical variables. By using 
the Hamilton equation of motion we can infer that the multipliers and are indeed equal 
to eo^ and They are dynamically undetermined pure gauge multipliers. 

In the PGT we consider a nine-parameter Lagrangian which is of the R + T^ + R^ type. It 
is well-known that the “most dynamical” case of the theory, when the parameter combinations 
are not critical, is of no physical importance due to the presence of ghosts and tachyonsB, one 
will expect there are “primary if-constraintR {P\G). Sec. 3.2 is devoted to the subject. Here 
we outline the procedure as follows: The hrst step is to hnd all the values of parameters which 
diminish the rank of the Hessian matrices d‘^C/dca^deiR and d’^C/dRa^’^dFR^ — the critical 
values of the parameters. Such values of parameters result in some torsion modes being frozen. 
We will show all of the possible primary constraints, which appear when the parameters take 
on the critical values. They will be written in the form of “if”-constraints, which automatically 
drop out from the theory when the corresponding critical values of the parameters are not 
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fulfilled. Finally, we will find the expression for the super-Hamiltonian which is valid for 
all values of the parameters. 

3.1 Total Hamiltonian 

Before we proceed to obtain the explicit form of the canonical Hamiltonian density and (pA, it 
is useful to define the decomposition of the related variables and functions. The components 
of the unit normal n to the =constant hypersurface, with respect to the orthonormal frame, 
are given by 

-e° 

A vector, e.g., can be decomposed into the orthogonal and parallel components with respect 
to the orthonormal frame indices: 




= -Vi_n^ 14 , 

(18) 

14 

III 

(19) 


= + n^rF). 

(20) 


One can extend the decomposition to any tensors with orthonormal frame indices. The lapse 
and shift functions can be written as 


N 

iV“ 




00 


= -n.eo^, 


9 


Oa 


^00 ^0 , 


( 21 ) 

( 22 ) 


and e = —NJ, where (—J) is the determinant of the 3-metric. 

To construct the canonical Hamiltonian density, we decompose the torsion and the curvature 
tensor as 


T — 

■L u, 


= 2n[rR^±^’' + RrU^'', 


(23) 

(24) 


so that and Rto^’' are independent of velocities. Defining the convenient “parallel” canon¬ 
ical momenta 

^ cr^a ^ (J /oc:\ 

^ IJ. — ^ ^1 ^ iMU — ^ /lui \2^) 

which satisfy ti"" nUa = 0, nyn„ = 0, we obtain 


fii/' 

deC 


(1 — Sq 


= J- 


dC 


deC 


The canonical Hamiltonian density. 


^ 


= J 


dC 




H = 7r“ e -t- -7r“ F - eH 
' ^can ~ 2 ^ ^ ^*-'7 


can be rewritten in the so-called 


Dirac-ADM fornilH 


(26) 


(27) 


Hc„ = ivHi + 1V“W„ + + dj>\ 


(28) 










7 


which is linear in N and iV“. The other quantities are given 


n± 

= + JC- 

(29) 

Ha 

= + \^\uRar - ea'^VfevrV, 

(30) 

Rfiu 


(31) 

pa 


(32) 


With the help of Eq. (p6D we can rewrite the expression in the parentheses in Eq. (p9D as a 
function of and Thus it is independent of the unphysical variables and 

represents the only dynamical part of the canonical Hamiltonian density. By combining ( 0 ) 
and (1^) the total Hamiltonian is 

Ut = NH^ + N’^Ha + ^ro^'^n^u + daV^ 

+ u^cl)A. (33) 

According to Dirac’s general arguments, there have to be (at least) ten hrst-class constraints 
in the PGT, as it is invariant under the ten-parameter local Poincare group. By utilizing the 
total Hamiltonian density (^^ in the consistency conditions for the sure primary constraints 
(pA|) we obtain the sure secondary constraints (SC) 


n± 0 , 0 , 


T-{ PC 


0 . 


(34) 


7-fx and Tia represent the generators of the orthogonal and parallel =const hypersurface 
deformationsi^@’0. And represents the generators of the local Lorentz transformationSil’ii. 
They are all hrst-class constraints, so that their consistency conditions are trivially satishedS. 
Only the super-Hamiltonian Tij, is involved in the dynamical evolution, the super-momenta TYq 
and Lorentz rotation parts are kinematic generators, consequently we concentrate on 
when consistency conditions are calculated. 


3.2 Primary if-constraints and super-Hamiltonian 

It is necessary to understand the relation between the canonical momenta and the velocities 
before one can make 7 Y_l more apparent. The torsion momenta can be decomposed into 
four algebraically irreducible parts: 

T T TTjjjj -|- ~rij[y7i, (35) 

where the notations refer respectively to the antisymmetric, symmetric-traceless and the trace 
parts. Manipulating the dehnition of the torsion momenta (0), the following relations between 
the different parts of the canonical momenta and the corresponding parts of the velocities T^-py 
are found: 

+-^{<^1 — o.2)Tp =-^{2ai + a2)T±pA, ( 36 ) 





IIV 


fj,l/ 

T 


+ ~ o.3)Tjiu± — -^(“i + 2a3)Tj_[jjj7], 


rp 

~77r -L I lTiT/\ 


J 


P 




j p 


(37) 

(38) 

(39) 


where T-p = Tpjj^, and a tensor with two indices contained in the bracket () denotes that the 
tensor is symmetric-traceless with respect to the two indices. If the parameters take on any 
of the critical values: 2ai -|- 02 = 0, ai -|- 203 = 0, ai = 0, or 02 = 0, the Lagrangian becomes 
singular with respect to some velocities (or equivalently, with respect to thus one 

A 

obtains the following PIC: (f)p± ~ 0, ~ 0, 0^^ ~ 0, and/or </> Ri 0, respectively. 


The curvature momenta can be decomposed into six irreducible parts: 


'^apLU 


and 




(40) 


flul. - 




TT: 


<7^U 


g ^apivl. 


^ [7i^]c 




a[iiv] 


Identifying the irreducible parts of the curvature momenta ( 0 ), one hnds 

= — - (&2 ~ = —(&2 + b ^ yR ^ o , 

J K K 


— '' M 


ayLV 


H —{bi — byRp± — - (&4 + b^Ripz, 

J Hi Hi 

+ -(61 ' 


J 


(t>i. = 

A 

yjiu± = 




Jl — b2)Tija7W± — —{bl + b2)^±wjiu^ 
K ' K 

^ ^ - 7^(64 - b,)R = \h, + 

J P 2k k 

A 

~r 

'^Jw± 

~r 


(&2 — byR ^ pp ] — —(62 + &5)-R±[77i7]±, 


K 


K 


— {bi — byRfjij^'^ — —{bi + &4)-R±(7n7)±, 
Hi Hi 


where ^Ro± := pR_lo := e^"''’'^Rrpua, Rptr ■= Rwjut", and R:= Rp 

gument as used above, for various degenerate parameter combinations one can obtain any of the 


(41) 

(42) 

(43) 

(44) 

(45) 

(46) 

(47) 

(48) 

By a similar ar- 


six expressions of (|43| - ^8D as PIC’s. The relations between the critical parameter combinations 
and the constraints are summarized in Table |l|. 

Blagojevic and Nikolic found that, for the generic PGT, there are only primary and sec¬ 
ondary if-constraints (no tertiary constraints). The Poisson bracket (PB) between either (i) 
appropriately paired primaries, or (ii) a primary and the secondary generated by its preser¬ 
vation, was found to be generally non-vanishing, hence they form a second-class pair. More 
specihcally, the value of these PB’s are just the (generally non-vanishing) constant “mass” pa¬ 
rameter combinations of Table |I| plus some held-dependent terms of nonlinear origin. In Sec. 5 
we will discuss the effects on the constraint classihcation due to these nonlinear terms. 
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Table 1. Primary if-constraints, critical parameter values and masses 


JP 

Kinetic Parameter 

Combinations 

Constraints 

Mass Parameter 
Combinations 

0+ 

(a) 

(b) 

02 

64 + bQ 

0, 

0-L, 

X 

x± 

oQ) 2 ao - 

1 - 02 


(a) 

di + 2a3 

077F> 

A 

X-flU 



1+ 

(b) 

1>2 + ^5 

A 

077F_L! 

A 

Xth/J- 

Ol — Oo, 

^0 , 
y +03 

2+ 

(a) 

Ol 

07717: 

Xjli7 



(b) 

bi + 64 

077!7_L: 

X-jiui. 

00,01 — 

Oo 


(a) 

2ai + 02 

077-L) 

Xtz± 



1" 

(b) 

64 + 65 

la. 

1 

Xj! 

Ol — Oo, 

2oo H“ 02 

0- 


62 + bs 

"0, 

A 

Oo , 

y +03 


2" 


bi + 62 

"0^7777: 


Ol — 02 



In order to treat all such possibilities in a concise way, the singular function 


X 


X 

0, X = 0, 


(49) 


was introduced. The PIC’s in the total Hamiltonian density (|33|) can be given in the form 

u^(t)A = {u ■ 0)'^ + {u ■ 0)^, (50) 

where 

{u ■ 0)^ = [1 - A(2ai + a 2 )]u'"-^(j>jzi_ + [1 - 

+ [1 — A(ai + 2a3)]M^^0jjj7 + -[1 — A(a2)]M0 (51) 

and 

(m ■ 0)-^ = -[1 — A (62 + &3)]ht'0 + -[1 — A(6i + 

+ [1 — A (64 + &5)]m^0^ + 2[1 — X{b2 + ^5)]^^^^07II7± 

+2[1 — A(6i + &4)]'tt^^^0jjj7^ + -[1 — A (64 + bQ)]u'^(l)±. (52) 


The super-Hamiltonian Hi_ in Tican (1^) then turns out to be of the form 

nj _ H-ZT I 'XjR 

/lJ_ — ^_L ' ^J_5 


( 53 ) 
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with 


Hi 


/ 

2 


3A(2cii “h CI2) 
2ai + Qj 2 


m = Jk 


Ql 

A(&2 + h) p,2 


czi + 2a3 

+ JC^ - n^^VaTT^ 


_ 24(62 + 63) 

A(6i + 62) T 

3(61 + 62) 

A (61 + 64 


+ 


+ 

A(Q2) ,2 

3(32 

A(64 + 65) 

4(64 + 65) 


/i5 


(54) 




2(61 + 64) 




2(62 + 65 
A(64 + 60) 
6(64 + be) 


+ JC 


R 


(55) 


where 


= 


12/2 L 


(2ai + 

+2(ai — — 2(ai — a2)T-jxT^ 


-c^Rr^R^ - C3(E^E^ + Rji^R^^) 

-c,R^R^ - cstf - + A, 


2P- 


(56) 


(57) 


here A is the cosmological constant. The relations between the constants q’s and the 6j’s are 
given by 


Co = 

24^ (2^1 + 362 + 63), 

Cl = 

CO 

1 

^o 

^ 1 ^ 

1 

C2 = 

24^(261 362 + 63), 

C3 = 

T( 5 i + *2 - 64 - h ), 

C4 = 

T( 54 - 52-54 + W. 

C5 = 

( 26 i 364 + 60). 

24 k 


(58) 

We now can use the formulations above to construct the total Hamiltonian under different 
critical parameter values. For example, if 2ai+a2 vanishes, this results in the relation in Eq. (^) 

^ 0. Therefore [1 — A(2ai + a 2 )]vP'^(j)-p± = vR-^cjyjzR 


becomes the constraint (j>p± = —^ + yyFj 

tJ L 


in the torsion if-constraint part (pTD; at the same time, —^= 0 in 
by referring to the singular function (|49|). On the other hand, if 2ai + 02 7^ 0, the relation (^) 
stands and ^ xp-i- _ -^- 4>T[±(jR^ and [1 — A(2ai + a 2 )]u^'^(j)-n± = 0. This 

2ai + 02 2ai + 02 

process can be performed term by term with different critical parameter combinations. It is 
very benehcial for the analysis of the PGT with general or specific values of the parameters. 




























4 The Linearized Theories of PGT 
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It is difficult to tell whether the PGT is viable under physically reasonable circumstances 
because of its complications (Indeed deep analysis has still left the issue obscure). Accordingly, 
the weak field approximation approach was the first and preliminary method applied to the 
PGT to study its properties. Linearization gives us insight into the viability of the propagating 
modes in the PGT. In the linearized theory we usually impose the physical requirements that the 
propagating modes carry positive energy and do not propagate faster than light. Theories with 
higher spin interactions are notorious for having such problems^. It was noted that interactions 
involving spins > 1 were especially vulnerable to having their characteristics outside the metric 
light cone, thereby permitting acausal modes of signal propagationlii. But linearization gets rid 
of nonlinear terms. One can tentatively forget the complex and focus on the structure of the 
theory. It is efficient for quickly eliminating the obviously unht theories. 

In the linearized PGT, in addition to the graviton there are possible propagating modes 
with spinP^'^^^5^ = = 2'^, 2~, I"*", 1“, O’*", 0“. The requirements for the possible propagating 

“massive” modes are essentially v < c (no tachyons) and E > 0 (no ghosts). Sezgin and van 
Nieuwenhuizenl’il studied the linearized PGT in the early 80’s. They applied the so-called 
“spin-projection operator formalism”! to the PGT and discovered that at least three of the 
massive modes have to be eliminated. By their analysis there are twelve such solutions. Their 
results are listed in Table |^. Hayashi and Shirafuji used the weak held approximation and 
the Klein-Gordon equation to investigate the particle spectrum of PGt!. They gave tables 
of all the maximal dynamic torsion modes. Their analysis conhrmed that these requirements 
preclude the possibility of more than 3 dynamic propagating torsion modes. Gonsequently, 
there are many distinct cases with the torsion modes showing different behavior: some modes 
frozen and some dynamic. 


Table 2 . Parameter choices for ghost- and tachyon-free gravity Lagrangian deter¬ 
mined by Sezgin and van Nieuwenhuizen. 



Parameter choices 




Particle content 

(1) 

ai -1- 203 

= bi + b 2 

= 64 

4 - 65 

= 0 

2+,0+,0- 

(2) 

Ol -|- 203 

= bi + b 2 

= 61 

+ 64 

= 0 

i-,o+,o- 

(3) 

Oi -|- 2 o3 

= 6i -f 62 

= 64 

4 - be 

= 0 

2+,l-,0- 

( 4 ) 

2ai + 0-2 

= 61-1-62 

= 61 

+ 64 

= 0 

i+,o+,o- 

( 5 ) 

2 ai + 02 

= 61-1-62 

= 64 

4 - be 

= 0 

2+,l+,0- 

(6) 

di — bi 

- 62 — 64 -I- 65 - 

= 0 


i+,o+,o- 

(7) 

ai = bi 

- 62 = 64 -|- 63 = 

= 0 


i+,i-,o- 

(8) 

di + 2 a 3 

= 62-1-63 

= 61 

+ 64 

= 0 

2-,l-,0+ 

( 9 ) 

2ai + d 2 

= 62-1-63 

= 61 

+ 64 

= 0 

2-,l+,0+ 

(10) 

2 di + d 2 

= 62-1-65 

= 61 

+ 64 

= 0 

2-,0+,0- 

(11) 

Ol = 02 = 

= 61 -1- 62 = 

= 0 



i+,o- 

(12) 

bi +b 2 = 

: 64 -|- 65 = 

= 62 4 

- ^5 = 

0 

o+,o- 


The linearized method was also applied to the cases with “massless” modesB0S. The 


propagating massless modes can then have gauge freedom. Kuhfuss and Nitsch! obtained 


some 
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very restrictive PGT parameter restrictions by using the spin-projection operator formalism 
and a “no fourth order pole requirement”. Their restrictions are 


A — 2(ii -|- 0,2 — T 2o2 — 0, 

+ &4 = &4 + &6 ~ ^2 T ^5 = 0, (59) 

Oi = Oo > 0, 62, ^3 G R (arbitrary). 


This parameter choice basically activates the spin-0“ and spin-2“ modes (refer to Table |l|) and 
suppresses all the other modes? These two propagating modes are both massless. 

The idea of “second-class pairs” of constraints introduced in Sec. 3 is useful in analyzing 
the linearized PGT. It shows that the Hamiltonian analysis is in complete accord with the 
linearized theory analysis of the propagating mode^SS. The cases where one or more of 
the “masses” vanish have also been analyzed by Blagojevic’s group using the Hamiltonian 
procedure—although only to linear orderly. The results are in accord with the linearized theory 
analysis of propagating massless modes and their associated extra gauge symmetrieJiBilS. 

It is satisfying that all of the aforementioned results are essentially consistent with one 
another, however, we should beware of the limitations of linearization. Although we may obtain 
a relatively quick understanding of the system via the linearization of the theory, the result from 
the full theory, especially a highly nonlinear system, does not necessarily keep the same traits 
as its linearization. This fact has made us suspicious regarding these pioneers’ proclamations. 
Especially, we view the conclusions of linearized theory for the massless modes as rather likely 
requiring qualitative revision from nonlinear effects and hence we regard them as much less 
hrm than those for the massive modes. Gonsequently, we don’t know how much conhdence 
to place in their conclusions. Our reservations concerning the linearized propagating massless 
modes in no way affect the main conclusions from the survey of the linearized PGT. However, 
it seems likely that including the nonlinear effect terms would modify the Poisson brackets of 
the constraints, converting many of the vanishing brackets to non-vanishing ones. In this way 
a constraint which was hrst-class, with vanishing (on shell) brackets with all other constraints, 
would move to the second-class indicating that a gauge symmetry of the linearized theory was 
broken in the nonlinear theory. Hence we expect that the present understanding of “massless” 
PGT modes and their extra gauge symmetries, being based on a linearized approximation which 
appears inadequate, will require substantial revisions. 

Aside from nonlinear effects, the analyses of the linearized theories of PGT offers the identi- 
hcation of the torsion modes, the critical kinetic parameter combinations, the mass terms and 
other useful relations. The critical kinetic and mass parameter combinations are given in Table 
1^. The necessary signature conditions for no ghosts and no tachyons are presented in Table 


5 Hamiltonian Analysis of Propagating Torsion Modes 

The works discussed in the last section showed that there are still many viable linearized theories 
of PGT. In the other words, the linearized theories with these parameter choices have good 
behavior in the weak held approximation. The next step is to consider if the full nonlinear 
theories with the same parameter choices still survive. It is quite important that the whole 
structure of the theory stays the same before and after linearization. (Theories that do not 


*Our understanding of the propagating modes from this parameter choice is different from Ref. 1 , further study 
is needed to clarify this point. 
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Table 3 . No ghosts and no tachyons conditions for the torsion modes 


JP 

Energy> 0 and Mass^ > 0 

0+ 

^4 ^6 ^ O7 

0002(02 + 2oo) < 0 

0- 

62 -t 63 < 0 , 

as + ao/2 < 0 

1 + 

^2 + ^5 > 0, 

(oi -1- 2 o 3 )(ai - ao)(a3 + oo/2) > 0 

1 - 

G + ^5 < Oj 

( 2 ai -1- 02 )(ai — Oo)(a 2 + 2 oo) < 0 

2 + 

Oo < 0 , bi 

-1- 64 > 0 , ooOi(oi — Oo) > 0 

2 " 

bi+b 2 < 0 , 

oi — Oo < 0 


have this character are inherently nonlinear, and are much more difficult to analyze. We are 
very reluctant to seriously entertain such theories unless we have no alternative.) In particular 
we should have no change in the number and type of constraints, for that is what determines 
the number of gauge freedoms and the number of dynamic physical degrees of freedom. If the 
answer is affirmative, the theory can be smoothly transformed from nonlinearity to linearity 
and vice versa with its properties unchanged. If the answer is negative, the properties of the 
theory would change fundamentally upon linearization. In fact it has been argued that such 
changes in the properties are accompanied by the occurrence of serious physical problems, 
in particular anomalous characteristics@’0. Hence a sound case of the PGT with “good” 
parameter choices should not exhibit such a change of properties behavior when passing to 
the weak held approximation. Thus this becomes a criterion that determines if some specihc 
parameter choice is good enough, or whether it should be regarded as non-viable. 

We are going to investigate the viability of the individual simple propagating spin modes 
and two better-behaved cases. Here “simple” indicates that all parameters will be set to vanish 
except the key ones. Such severe conditions could be generalized. For example, in order to 
suppress all except the O’*" and 0“ modes, according Table 0. our condition could be relaxed to 
require only h\ = —62 = —&4 = h instead of our severe “simple” choice 61 = 62 = = 0. 

There are two reasons not to adopt the “relaxed” versions: (a) it would much complicate the 
analysis, (b) the relaxed parameter choices do not essentially enhance the viability of the theory; 
instead, they could even ruin a case which is well-behaved under the simple version. We can 
readily see this from the primary Poisson matrices of all the cases. 

It is not the purpose of this paper to examine viability in the full nonlinear domain of all the 
“viable” linearized theories of PGT. It seems appropriate to hrst try to illustrate the influence of 
nonlinearity in these theories and to investigate the effects in order to get some understanding. 
Nevertheless, we believe that the phenomena resulting from nonlinearity in the cases discussed 
in this paper are representative for all the other cases. We would like to emphasize the effects 
of nonlinearity on the change of the structure and gauges of the theory, so we will focus on the 
change of the type and/or number of the constraints and the consistency of the whole system 
caused by nonlinearity. These observations will be solid enough for an understanding of the 
differences before and after linearization. 

In our investigations we found that the “constraint bifurcation” phenomena is a key evidence 
which determines if there is an anomalous characteristic. Whether the Poisson matrix has 
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constant rank or not, independent of the field valnes, is the criteria for knowing whether the 
phenomenon of constraint bifnrcation occurs. A held dependent rank can happen only if the 
constraints are nonlinear. Strange behavior may occur as one approaches a point in the phase 
space where the rank of the Poisson matrix changes. At such points the number and type of 
constraints (and thus the number of gauges and physical degrees of freedom) are not constant. 
This phenomenon signals tl^ presence of an acausal propagating mode in the system and 
is unacceptable theoretically^. Whether the Poisson matrix is singular can be judged by the 
positivity of its determinant. Since the matrix is usually big the concise form of its detominant 
is hard to obtain, so we appeal to the aid of the symbolic computing software REDUCE^ for the 
calculations. From the result associated with all the constraints we can decide the positivity of 
the determinant and thus the possibility of constraint bifurcation. 

5.1 Simple spin-one and spin-two cases 

The method of the analysis is simple conceptually: list all of the constraints and compare 
the differences of the structure before and after linearization. The Poisson matrix is formed 
containing the non-zero Poisson Brackets (PB), then checked for change under linearization. 
Due to the complications, the secondary constraints (SC) will not be expressed explicitly. 
Instead we will reason out the differences based on the result of the SC’s linearization in 
certain situations. 

5.1.1 The simple spin-l+ case 

^]iu± corresponds to the spin-l’*' mode. The linearized theory indicates that the theory with 
only a spin-l+ mode can be obtained with the simple parameter choice: 


do 7^ 0, ds 7^ 0) ^5 7^ 0, 

di = d2 = 0, (60) 

bi = b2 = bs = bi = be = 0 . 

With such a parameter choice there are twenty-one PIC’s. The PIC’s are as follows: 


III III 

^1 d’^i 

p o 

III III 

k b 

-^-0, 

^_ 

7" ^ 0, 


(61) 

TT 3do „ 

= 

J 



The non-zero PB’s are 






o A 

2 



(62) 


^xx' A 

2 J2 


(63) 


cfbxx' pal. 
^ J2 A 

A 


(64) 

{0±,0'} = 

^xx' 6(Zo 



(65) 

J ’ 
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.P 

(66) 

r<PaW^(l)¥} = 

T7 ^ \ 

j2 \Vrl7'^'aJI±) 

(67) 

{>57717,0^} = 

1 

2 P 

'^JlTlv{p'^T)a o^Jl{p‘^f)v) 

. ( 68 ) 

{4>-pu±j 4>w±} = 

1 ^ 

2 J2 

(69) 

{0JZI7±) 075^} 

1 A ^0 

J pPiiJ.'nvya ■ 

(70) 


With these PB’s we can construct the primary Poisson matrix with the following form: 


13 5 115 5 







0 

07I-L 

(^w 

0 ± 


0m;7± 


1 

{ 

0 

0 

0 

0 

h 

0 

0 

0 

3 

{ 

0 jr± 

0 

71 

0 

71 

TT 

71 

TT 

5 

{ 

iw 

0 

0 

71 

0 

0 

r]+7r 

TT 

1 

{ 

0 ± 

V 

TT 

0 





1 

{ 


0 

TT 

0 



0 


5 

{ 

077I7± 

0 

71 

r^+vr 





5 

{ 


0 

TT 

71 






(71) 


where r^’s represent the constant terms and tt’s represent the terms involving the variables. The 
details of these terms can be found in Eq. O-O- The super-Hamiltonian density is deduced 
by applying the parameter choice (^) to (^). It is 



hT* 


+ nr, 


(72) 


where 


n 


Rl+ 


n 


T 1 + 


1^ 

2 


1a 


K- Tl-f} 

265 J 




A 

'^jn7± 


265 J 




Oo 




3/2 

2ci3 J 


I rj^fluA- 


Q 3 


7^2 - n^V„7r^ 


/!• 


(73) 

(74) 


Within the phase surface of primary constraints, the super-momenta (PP]) and the Lorentz 
rotation parts (|3l|) lead to the following constraints: 




n 


1 + 






+ R, 


^ua± 


Ifii/ a J I " J 

1 ^ rjiV ^ Tfi _ ^ ^ ^ 

A_ ^ ^ .. Aj7_ 

^0 rji I ^ 1 / r (7V7 ^ ^ 

J ^ + j7i[ji^u]n - ^u-j- 




ao 

¥ 


R 


■)i± 


0 , 


0 , 


(75) 

(76) 
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QQ rp I ^ ^ T" 


J 


2 / 2 ' 


1 rcr V7 ^ 
"2^ j 


S3 0. 


(77) 


It is essential to compare the nonlinear spin-l+ case with its linearized version. At hrst let 
us describe briefly the analysis of the linear spin-l’*' case. After linearization, the remaining 
non-zero PB’s turn out to be Eqs. (|65D , ([70|) in which only the zero-order (constant) terms 
exist, i.e.. 




(78) 

(79) 


where the superscript b indicates linearization. With the linearized PB’s we can construct the 
constant primary Poisson matrix: 


13 5 11 5 


5 








aflu 


1 { 

3 { 4x 
5 { 

1 { A 
1 { 

5 { 

5 { 


o 

-V O O 
0 0 0 

O O -V 
0 0 0 


h o o 

0 0 0 

O O h 

o 


o 

o 

o 


(80) 


where rj represents the constants in Eq. (ffSl) and (^). The consistency conditions guarantee 
that the linearized Lagrange multipliers, i.e., u^, u^-^, and are obtained. The nine 

linearized SC, also show up when the consistency conditions of the linearized 

constraints 0^_|_, 'fb'’ and are required. Furthermore, the linearized SC x^±: and 

will become the counterparts to the second-class pairs of 0 ^_|_, and respectively accord¬ 
ing to Ref. 23. These second-class pairs determine the multipliers hi*’ and and the algo¬ 
rithm of the linearized case is terminated. In the PGT, there are forty dynamic variables coming 
from the sixteen tetrad components and the twenty-four connection components. The number 
of total variables counts eighty because the same number of canonical momenta accompany the 
variables. Nonetheless, constraints eliminate many unphysical variables. It is known that the 
super-Hamiltonian super-momenta Tia, and the Lorentz generators are ten first-class 
constraints. There are also ten sure first-class primary constraints. The total twenty first-class 
constraints (because of their gauge nature) offset forty variables. The degrees of freedom of the 
linearized case can be counted as |(80[total] —40[gauge] — 21 [PIC] —9[SC]) = 5 = 3[1+] -|-2[GR], 
a massive spin-l+ field in addition of the usual graviton. For the sake of “ghost-free” the posi¬ 
tivity of the kinetic energy density demands 03 < 0 and b^/n > 0. 

Notwithstanding the good behavior of the linearized case, the nonlinearity of the full spin- 
l’*' theory ruins its whole legitimacy. Generically, the Lagrange multipliers u, «■*■, 
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can not be determined alone becanse the nonlinear terms are engaged (see Eq. O)- Here we 
assnme tentatively that there is no inconsistency appearing in the system dnring the process. 

The consistency conditions of the connection PIC’s, 0 _l, and ^ 5 ^, will determine the 

multipliers of the tetrad PIC’s, u, and u^’^, and give three SC’s, X[ 3 ]. Principally, the 
consistency conditions of the tetrad PIC’s and the SC’s ypj will determine the multipliers, «■*“, 
'ti, and according to the experience of the linearized case and terminate the process. 
The degrees of freedom of the full nonlinear case count as |(80 — 40 — 21 [PIC] — 3[SC]) = 8 = 
2 + 3[1’''] + 3[1“]. This shows that the massive spin-l“ held is excited as well as the massive 
spin-l"*" held since the other modes are already suppressed. The result is unacceptable from the 
point of view of the positivity of the kinetic energy density in the Hamiltonian density (refer 
to Eq. (P). 

One can learn the major diherence just from the primary Poisson matrix. The rank of the 
linearized primary Poisson matrix (^) in this case is twelve, a constant, whereas the rank of 
the full nonlinear matrix ([7TD varies. Its rank is eighteen generically but it can be less than 
eighteen depending on some specihc numerical values of the variables which could degenerate 
the system. We verihed this point by using the symbolic computation software REDUCE to 
calculate the determinant of the matrix The result shows that the determinant is not 

positive dehnite, which indicates that whether the system is singular depends on the numerical 
values of the variables. This situation shows the phenomenon of constraint bifurcation. It could 
even happen as the system dynamically evolves. Consider all of the possible initial values that 
give generic rank to the Poisson matrix. We see no mechanism that would guarantee that every 
one of these cases could not evolve to held values where the rank changes. Because of these 
features we regard this parameter choice as physically unacceptable. 

5.1.2 The simple spin-l“ case 

We next examine the simple spin-l“ case and hnd that it has similar problems as in the simple 
spin-l'*' case, tt^ corresponds to the spin-l“ mode. The linearized theory indicates that the 
theory with only the spin-l“ mode is obtained if only oq, 02 and 65 do not vanish. Under this 

A 

parameter choice 0^ instead of 0 and 0^^ become PIC’s and the other PIC’s remain just as 
for the spin-l’*' case. The non-zero PB’s are 



J2 

(81) 

{ 0 ±, P'-JUt] 

^xx' /\ 

J2 

(82) 

4>'tp} 


(83) 

~ A 

^xx' A 
~ J2 

(84) 

0 ^} 

= rhs of Eq. ( 68 ), 

(85) 


= rhs of Eq. ([70|) . 

( 86 ) 
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By utilizing the PB’s above the primary Poisson matrix has the form 

1 3 3 5 1 1 5 5 


1 

3 

3 

5 

1 

1 

5 

5 


{ 


A 

fw 

iw 

<P± 


0 0Tr± 0; 


V 

O 

O 

O 


The super-Hamiltonian density is 
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TT 

71 

71 
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71 

71 
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0/i!^± 


V 

0 

0 

0 

71 

71 

71 

71 

71 
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71 

71 

0 

0 

7 ] + 7 T 

71 


= n 


O 

O 

?7 + vr 
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o 


(87) 


( 88 ) 


where 


'XjTl _ 

^_L 111. 


3/2 7l^^^ 

2,0,0 J 


+ T^ 


/2 TT^ 

+ V-V - n^V„7r' 


6a2 J 




(89) 


The constraints derived from the super-momenta and the Lorentz rotation parts are: 




1 

Rjn7a± '^^RjTu'^ 


Oo 7l’'± 

^2 “1“ j 


iv-— ^ 0 

3 


(90) 


n 


1+ 

n± 


IJ,U± rpT7 _ _ ^0 rp 

J J P ^ 


2 _ 

+ jl[-pVu]n'^ - 5/Vu^ ~ 0 , 


R-w 


On ^ 1 ^ zi; 


-^0 [7lVx7]— - JTT [JZ|_ 


VT^Uh 


0 . 


(91) 


(92) 


By using same argument as in the simple spin-l’*' case, it is obvious that nonlinearity remains 
untamed. Linearly the degrees of freedom are |(80 —40 —20[PIC] —10[SC]) = 5 = 3[1“]-|-2[GR]. 
The spin-l“ held is obtained as we expected. But nonlinearly it becomes ^(80 — 40 — 20[PIC] — 
4[SC]) = 8 = 3[1“] -|- 3[1+] -|- 2 . The appearance of spin-l+ allows the possibility of negative 
kinetic energy and hinders the success of the case. The possible constraint bifurcation is noted 
and its consequence is similar to the spin-l+ case discnssed above. 

In the aforementioned spin-one cases onr analysis shows that neither the simple spin-l+ 
mode nor the simple spin-l“ mode can be separated from each other. To see this more clearly, 
we can consider a wilder parameter arrangement: Oq and 65 non-vanishing only. This choice 
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has the consequence that all tetrad variables are not dynamical; they turn out to be the tetrad 

PIC’s: (j), 0JJJ7 and The other PIC’s remain the same: ^(p, (p±, and (pjiu±- 

We hnd that generically the degrees of freedom count as |(80 — 40 — 24[PIC]) = 8 = 
3[1“] + S)!"*"] + 2. The only difference between this “simpler” case and the simple spin-one 
cases discussed above is that the SC’s induced from the PIC’s in the simple spin-one cases 
have essentially transformed into PIC’s coming from the parameter choice in this “simpler” 
case. However, the spin-one modes still stick to each other to crash the theory because of the 
negative kinetic energy. 


5.1.3 The simple spin-2 case 


With the understanding of the simple spin-one cases the spin-2’'' mode which stands with the 
parameter choice of only oq, oi and bi non-vanishing is not expected to have good physical 
propagation because it is entangled in a more complicated nonlinearity. Instead we examine 
the simple spin-2“ case. Although it eventually becomes futile, it is illuminating for our next 
case. According to Table |I| corresponds to the spin-2“ mode. We make the specihc 
parameter choices: 

0-0 7 ^ 0 , bi ^ 0 , 

ai = 02 = as = 0, (93) 

62 = 63 = 64 = 65 = fee = 0. 

The PIC’s are as follows: 





The non-zero PB’s are 
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and the primary Poisson matrix has the form 
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The corresponding super-Hamiltonian and the reduced constraints are 
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(108) 


Guided by our earlier experience, we readily realize that after its linearization the remaining 
non-zero PR’s turn out to be Eq. (P^, ^ , |10^) in which only the zero-order (constant) terms exist. 
Then the consistency conditions guarantee that the linearized Lagrange multipliers, i.e., 

.yb/ia, :^b)iu are obtained. The six PIC’s, and 0^, and the six SC’s, and 

X^, which are derived from the consistency conditions of ^0^ and 0^ are second-class pairs, 
respectively. These second-class pairs determine the multipliers and and the algorithm 
of the linearized case is terminated. The degrees of freedom of the linearized case are counted 
as 2(80— 40 — 20[PIC] —6[SC]) = 7 = 5[2“] -1-2, a massive spin-2“ held in addition to the usual 
graviton. 

On the other hand, the full spin-2“ case does not follow this track. Generically, «■*■, u^, 
can be determined. Since *0 does not commute with 0^^, the consistency condition of 
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(p-pjj can determine and give fonr secondary constraints xi^]- The consistency conditions of 

A ^ 

the connection primary constraints, (j)±, ^(j), (p-pjj and (p-p, leave fonr of the tetrad Lagrange 

mnltipliers, u, and nndetermined. From the experience of the linearized case, we 

infer that the consistency condition of the fonr xi 4 ] will allow the fonr nndetermined tetrad 
Lagrange mnltipliers to be worked ont and thereby terminate the process. The degrees of 
freedom wonld then count as |(80 — 40 — 20[PIC] — 4[SC]) = 8 = 5[2“] + 1 + 2. It appears that 
one degree of freedom of the spin-2+ mode is excited as well as the degrees of freedom that 
we expected. This is rather suspicious; a more careful analysis should be done to clarify these 
details. However, the main point is certain: once again the nonlinear effect makes the structure 
of the full simple spin-2“ case different from its linearized case. 

The ill-behaviors in the cases above show that the existence of a single higher-spin mode 
is virtually impossible. Generally, in these cases there are extra degrees of freedom and the 
phenomenon of constraint bifurcation (compared with their linearized theories). Here the extra 
degrees of freedom represent ghosts, i.e., propagating negative energy, and the phenomenon of 
constraint bifurcation leads to tachyons, i.e., faster-than-light propagation. Therefore, nonlin¬ 
earity makes the full theory qualitatively different from its linearization and the theory fails 
important theoretical tests. 

The constraints coming from the reductions of the super-momenta and the Lorentz rotation 
parts do not seem to be much help in determining the existence of the extra degrees of freedom 
and the constraint bifurcation. This is because these constraints can be used to eliminate 
some unphysical canonical variables instead of unphysical canonical momenta. But only the 
canonical momenta are involved in the phenomena in these simple cases. 

The reason we do not show the explicit forms of the SC’s in these cases is because they are 
very complicated and we can obtain a sufficient understanding with only their linearized forms. 
The detailed form of these SC’s is certainly needed in order to solve for the explicit form of the 
Lagrange multipliers. But the main point is that we know that this will be possible since they 
will not commute with all the primary constraints (in particular with their second-class pair 
counterparts). The complication of all the SC’s in the full nonlinear theories comes from the 
entanglements of the variables and the momenta. These complications will not decrease the 
possibility of constraint bifurcation for the same reasons that were stated in the last paragraph. 

5.2 Simple spin-0“ + spin-2“ cases 

Although the inconsistency of the simple spin-2“ case with its linearization renders the theory 
unsuccessful, the sign of the few unexpected degrees of freedom appearing after including the 
nonlinear terms is notable. It becomes straightforward to recognize the situation from Table |I[ 
In Table ^ only the PIC’s and do not have their “natural” second-class pair counterparts 
coming from the tetrad PIC’s; their second-class pair counterparts are the SC’s derived from 
their own consistency conditions, i.e., ^<j) and Thus the status of these two modes is 

unconcerned with the the tetrad parameter choice even in the full nonlinear case. This feature 
can be used advantageously; it will help in constructing some “more viable” simple cases. Here 
we show two parameter choices which give propagating ■+ and modes generically, albeit 
with drawbacks in the possible degeneracy. 
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5.2.1 The negative 62 case 

A 

According to Table |I| the parameter 62 appears in ^0(0“), and Here we 

make the tetrad parameters vanish in order to suppress the spin-l’*' mode. Therefore only the 
spin- 0 “ and spin- 2 “ modes propagate, the specihc parameter choice is 


7^ 0) &2 ^ 0, 

Ui = 0,2 — CL 3 — 0, (109) 

61 = 63 = 64 = 65 = 6g = 0. 


This choice leads to the PIC’s 
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The primary Poisson matrix is 
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Although this case seems more complicated, we found that the theory now has generically the 
same structure as its linearization. Briefly Eq. ( |111| , |113| , |118|) are the only non-vanishing PB’s 
after linearization because of their constant terms. Then the linearized multipliers 

and can be determined from the consistency conditions of the linearized 

PIC’s (j)^, 0^, 0^_|_, 0p 0^ and 0^_i_. Since 0^ commutes with all the other linearized PIC’s, 
its consistency condition gives the linearized SC’s, i.e.. 
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As expected appears in so and 0^ are a second-class pair. Therefore the SC’s 
can be used to determine the multiplier At the same time also suppresses the the 
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spin-l'*' mode and prevents the system from having the problem of a propagating negative 
energy mode. 

Generically the structure of the full theory is consistent with the one of the linearized 
theory despite the fact that the analysis is not as straightforward as the description in the 

A 

last paragraph. The key difference is that does not commute with the connection PIC’s 
nonlinearly. The consistency condition, 



'jLU' "^t} 


0 , 


(125) 


can give the SC XjuJi here the connection multipliers shown in Eq. (|125|) can be derived from 
the consistency conditions of the other tetrad PICs. Xjiu includes as the linear part and 


thus associates with 0^ to become a second-class pair. Consequently, the SC renders 0^_i_ 
non-dynamical. Now we have enough constraints to determine all the multipliers and hnish the 
process. The degrees of freedom are |(80 — 40 — 21 [PIC] — 3[SC]) = 8 = 1[0“] -|- 5[2“] -|- 2 in 
both the linear and nonlinear theory. The result indicates that the theory can be transformed 
into its weak held approximation without qualitative changes in the canonical structure. 

However the theory cannot entirely escape from the constraint bifurcation problem. By 
using REDUCE the possible vanishing of the determinant of the sub-Poisson matrix formed by 
the PB of 0^^, 0jjjj, 0JJ and 0^^ is affirmed, i.e., with certain numerical values of the canonical 
variables which satisfy all constraints the sub-Poisson matrix can be singular. This is not so 
good. Such an occurrence is thought to be a sign of the presence of a mode which propagates 
faster than hght0. 


5 . 2.2 The hi + case 

The other parameter choice available to activate only the spin-0“ and spin-2“ modes is to let 

A ~ 

all parameters vanish except Oq, hi and 63 . This choice gives the PIC’s 0, 0^j_, 0 jjj 7 , (j)±, 

A ^ 

0^_i_ and (j)p. The non-zero PB’s are 
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and the primary Poisson matrix is 


1 

- — hre 

.p [e 

^xx' T 
j2 ^WJIu- 


5 1 


0M± 0/717 0/7i7 0± 0/7i7± 0y 


(132) 

(133) 


TT TT 


1 { 0 h o o 

3 { 0y± ^ O tt r/+7r 

A U 

3 { 0/717 O h+TT TT 

5 { 0/717 TTTT- 

1 { 0± h O O TT 

A 

3 { 0^± O tt r/+7r tt O 

3 { 0/7 O h+TT TT TT 

The corresponding super-Hamiltonian and the reduced constraints are 
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Essentially the same argument as used in the negative 62 case can be applied to this case. 
The degrees of freedom are generically |(80 — 40 — 19[PIC] — 5[SC]) = 8 = 1[0“] -|- 5[2“] -|- 2 in 
both the linear and nonlinear theory. But once again it is readily seen that this case also almost 
certainly has constraint bifurcation and consequently is very vulnerable to acausal propagation. 

Two questions are raised by these results: (a) Could relaxed versions of the two spin-0“ 
-|- 2“ modes overcome the defect of constraint bifurcation? (b) Do there exist any parameter 
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choices which make the compound propagating modes nonlinearly viable? Concerning the 
first question, according to the explanation in the beginning of this chapter, more nonlinear 
effects will be involved in the relaxed versions. The Poisson matrix is supposed to be more 
complicated and even include the spatial derivatives of If will become harder to analyze the 
whole system and to convince people that it has the advantage of being free of any constraint 
bifurcation problems. To the second question, basically this is beyond our ability to give a 
definite answer. Maybe with some delicate tuning of all parameters there exist such nonlinearly 
viable cases of the PGT. (After all, unexpected miraculous cancellations have been found in 
other systems.) However, based on our experience with the PGT, we hold a very conservative 
attitude to the possibility of this happening. 

In order to convince ourselves that the cases dubbed “simple” in this section do show the 
obstacles that all the relaxed cases and many similar cases could meet in the Hamiltonian 
analysis, it seems worthwhile to investigate the most promising parameter choice coming, from 
the linearized PGT. In section 4 the parameter choice (|^ given by Kuhfuss and NitschlU 


was 


quite restrictive and supposed to be “viable”. Here we give a brief analysis of the theory with 

A 

this parameter choice. Referring to Table |I| we learn that the PIG’s in the theory are 

and 0^^. Under linearization one will expect that 0A and 0^_i_ are first-class due 

to ai = Oq. And 0^_|_, and 0^_i_ commute with one another. Subsequently x^_|_, x\. and 
Xjiu± are derived from their parent constraints. These SG’s and their parent constraints form 
second-class pairs and thus suppress 7rh_|_, and respectively. The degrees of freedom 
count as |(80 — 40 — 15[PIG] — 9[SG] — 6[gauges]) = 5 = 2[GR] -k 3. It is easy to discover that 
the three degrees of freedom come from the massless spin-0“ and spin-2“ propagating modes. 

Now we return to the full nonlinear case. In order to make the analysis simpler, we can 

A 

specify the parameters to be 6i = 64 = fee = 0. First, 0^^ and 0^_i_ will not commute with the 
connection PIG’s because of the nonlinear terms. Gonsequently we expect that the consistency 
conditions of all the PIG’s will determine twelve Lagrange multipliers and give three SG’s ypj. 
Using X[3], the remaining three multipliers can be obtained. The degrees of freedom would then 
count as |(80—40 —15[PIG] —3[SG]) = 11 = 2[GR]-|-l[0“]-|-3[l“]-|-5[2“]; all excited propagating 
modes have become “massive” formally. Nonlinear effects have excited extra dynamic degrees 
of freedom and destroyed the gauge freedom of the linearized theory. Moreover, if we require 
the kinetic energy density of the spin-2“ mode to be positive definite, it is inevitable that the 
spin-1 “ mode propagates with negative energy density. This unexpected mode with the wrong 
sign shows another defect of this parameter choice. 

Here we try to illustrate that the relaxed cases can be even more vulnerable to getting 
stuck in troubles. Essentially, the Kuhfuss & Nitsch’s case is in the same track as ours in 
Section 5.2. They concluded that only the spin-0“ and spin-2“ modes could propagate. We 
found from the nonlinear structure of the PGT that it is usually unlikely to suppress one single 
higher-spin (s > 0) mode without suppressing its counterpart if they form a second-class pair 
in the linearized theory. This realization forces the parameter choice to be the one which only 
activates the modes not belonging to any second-class pairs — the spin-0“ and spin-2“ modes. 

Maybe there exists still the thought that it might still work with different values of the 
nonzero parameters or in more relaxed versions. Until we do further more detailed investigations 
this possibility cannot be dismissed. Along this line one has to face more problems: (1) Because 
the connection PIG’s do not commute generally in the relaxed cases, more non-zero PB’s appear, 
which could excite extra degrees of freedom. Negative energy can be expected to appear in 
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most of those cases; (2) In general, the spatial derivatives of the 6xx' function could appear 
in the PB’s of the connection PIC’s. This greatly complicates the difficulty of analyzing the 
generic rank of the Poisson matrix. Even if these problems are settled down well there is the 
very possible phenomenon of constraint bifurcation with its attendant acausal modes waiting 
ahead! 


6 Discussion and Conclusion 

In this paper we examined the behavior of a restricted case of the MAG — the PGT, the 
local gauge theory of the Poincare group. People used to work on the linearized PGT and its 
initial value problems to clarify its viability, however they overlooked the signihcant influences 
of nonlinearity on the whole theory. In Sec. 3 we noted that two problems are produced by 
nonlinearity: (1) Constraint bifurcation-, the phenomenon is caused by the appearances of the 
nonlinear terms in the constraints which lead to a held dependence in the Poisson matrix. 
This can cause a bifurcation in the constraint chain. The number and type of constraints 
can depend on the hnld values. This phenomenon has been linked to acausal (tachyonic) 
propagation modes^^J. (2) Field activation: the phenomenon is that nonlinearity turns some 
original constraints in the linearized theories into held equations. This means that some held 
modes which should be frozen are activated. If the helds carry negative energy, it violates the 
“no-ghost” requirement. We regard these two phenomena as very important; our criteria is 
that they should be avoided. A good theory should keep the same dynamical structure before 
and after linearization; it should not be bothered by these two phenomena. We hnd it hard 
to imagine the consequences of a theory of a fundamental interaction which would exhibit a 
diherent dynamical structure as we passed from the strong held to weak held regime. 

In the last section we examined nonlinear ehects in the PGT which can switch the status 
of presumably viable parameter choices. The main tools are the Dirac-Bergmann algorithm 
of Hamiltonian analysis, and the if-constraint technique developed by Blagojevic and Nikolic 
(described in Sec. 3). The procedures we used are as follows: 

(a) Make the desired parameter choice according to the study of the linearized PGT. 

(b) Identify all the constraints in the selected case and calculate their Poisson brackets. 

(c) Glassify the constraints by the results of the Poisson brackets. From the consistency 
conditions hnd out all secondary constraints and the Lagrange multipliers. 

(d) Gount the degrees of freedom of the selected case. By knowing the corresponding lin¬ 
earized results, we can hgure out the meaning of each degree of freedom without further 
calculation. 

(e) Evaluate whether the Poisson matrix formed by the Poisson brackets of all the constraints 
can be singular at specihc values of the variables in order to determine any occurrence of 
constraint bifurcation. The computer algebraic software REDUCE can be applied to work 
out the lengthy calculations. 

Individual “simple” spin modes, which are well behaved and thus viable in the linearized 
PGT, were chosen to study their behavior under the full nonlinear considerations. For the 
spin zero modes, we previously found that they are essentially viable although nonlinear effects 
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obviously complicate the whole systems We noted that there is the possibility of constraint 
bifurcation, however, this can be avoided with specihc value choices of the parameters^. Here 
we concluded that the nonlinear terms simply devastate the viabilities of the spin one modes 
and spin two modes, since each expectedly frozen mode which has the same spin but opposite 
parity with the propagating mode in every case is activated by nonlinearity. Unfortunately, 
the nonlinear-activated mode propagates with negative energy which is physically unacceptable 
theoretically. 

There can be maximally three different propagating torsion modes in the linearized PGT 
as described in Sec. 4. Therefore we looked for the possibly viable multi-modes according to 
the classihcation in Ref. 7. Through many trial-and-errors only two “simple” cases of multi¬ 
modes which circumvent the difficulty of the activation of the unexpected modes are found, 
i.e., keep only the spin-0“ and spin-2“ modes alive and suppress all the other spin modes. 
Successfully the structure (and the degrees of freedom) of the full nonlinear theory remains the 
same as that of its linearization in these two case. This is a promising results. Even though 
we argued that the constraint bifurcation of the Poisson matrices could be expected to render 
them unattractive, we wish to note that that issue, and these modes more generally, are worthy 
of a deeper study. 

There are a few other special parameter choices meriting examiniation, e.g., those identihed 
by Katanaev which give solutions with vanishing curvature or torsion0. One that we specihcally 
considered here was Kuhfuss and Nitsch’s case, which was thought to be the most promising case 
for being viable. We found that the case changes its character drastically between the nonlinear 
and linearized versions. Both the phenomena of constraint bifurcation and held activation 
happen; they devastate the case. Considering this result, one could be pessimistic about hnding 
viable higher spin modes in the nonlinear PGT. But we don’t exclude the possibility, with a 
hne tuning of all the involved parameters. Once the case exists, the propagating modes are 
expected to be short-range, i.e., massive, just like the viable spin zero modes. 

It seems that only a few cases in the PGT can avoid the phenomenon of held activation. 
And almost all of these cases are plagued with the constraint bifurcation problem. This means 
that the nonlinear PGT is very likely to be qualitatively diherent from the linearized PGT in 
the number and type of constraints. The more general MAG theory includes the PGT as a 
subcase and thus already has these same problems. Moreover it seems extremely likely that 
nonmetricity itself will be vulnerable to the same kind of problems. 

Our analysis of non-linear constraint ehects in the PGT reveals an extreme difficulty in 
hnding a viable gauge alternative of gravity which had not been appreciated in the past. In the 
light of this, a further important consequence of our analysis is an enhanced understanding of 
the advantages and uniqueness of GR. 
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